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Abstract 


Mechanics for the analysis of laminated composite shells with piezoelectric actuators and sensors ^re 
presented. A new mixed-field laminate theory for piezoelectric shells is formulated in curvilinear 
coordinates which combines single-layer assumptions for the displacements and a layerwise 
representation for the electric potential. The resultant coupled governing equations for curvilinear 
piezoelectric laminates are described. Structural mechanics are subsequently developed and an 8 -node 
finite-element is formulated for the static and dynamic analysis of adaptive composite structures of 
general laminations containing piezoelectric layers. Evaluations of the method and comparisons with 
reported results are presented for laminated piezoelectric-composite plates, a closed cylindrical shell 
with a continuous piezoceramic layer and a laminated composite semi-circular cantilever shell with 
discrete cylindrical piezoelectric actuators and/or sensors. 
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1. Introduction 


# 

Smart composite laminates and adaptive composite structures with embedded piezoelectric sensors 
and actuators seem to combine some of the superior mechanical properties of composites with the 
additional capabilities to sense deformations and stress states and to adapt their response accordingly. 
Consequently, such novel and structures are receiving substantial research interest. Curvilinear shell 
structures are commonly used in aeronautical, aerospace, automotive and other engineering 
applications. Thus, adaptive curvilinear composite shells with embedded piezoelectric devices are 
among the structural configurations of most practical interest, yet, they are also among the more 
challenging to study both analytically and experimentally. The present paper presents theoretical 
foundations and mechanics for the analysis of curvilinear piezoelectric laminates and shell structures. 

Numerous theories and models have been proposed for the analysis of laminated composite beams 
and plates containing active and passive piezoelectric layers. Simplified approaches attempting to 
replicate the induced strains and electric fields generated by a piezoelectric layer under an external 
electric field have been proposed (see e.g. Lee and Moon (1989), Lee (1990), Lee and Moon (1990), 
Wang and Rogers (1991), Lazarus and Crawley (1989)). Variational methods and finite element 
models for piezoelectric solids have also been reported by Allik and Hughes (1970), Naillon et al. 
(1983), Tzou and Tseng (1990), and Ha et al. (1992). Layerwise theories were reported for infinite 
piezoelectric plates (Pauley (1974)) and for finite elastic laminated beams and plates with induced 
strain actuation (Robbins and Reddy (1991, 1993)). Coupled layerwise theories and finite elements 
for laminated composite beams and plates with piezoelectric actuators and sensors were also 
developed (Heyliger et al. (1994), Saravanos and Heyliger (1995), Saravanos and Heyliger (1996)) 
which consider the complete electromechanical response of smart piezoelectric plate structures under 
external mechanical or electrical loading. An exact piezoelasticity solution has been also reported by 

t 

Heyliger and Saravanos (1995) for piezoelectric laminated plates. 

Analytical formulations have been proposed for laminated piezoelectric shells. Dockmeci (1990) has 
reported theoretical work for the vibration of single-layered piezoelectric shells. Lammering (1991) 
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developed a Reissner-Middlin type shear deformable finite element for shells with surface bonded 
piezoelectric layers. Tzou and co-workers have presented approaches for the analysis of thin 
laminated shells based on Kirchoff-Love assumptions (Tzou and Garde (1989)), and theoretical 
formulations for shells which include Reissner-Middlin type of shear deformation theory and rotary 
inertia effects (Tzou and Zhong, 1993). Some of the previous approaches effectively utilize equivalent 
force/moment representations of the induced strain by attempting to replicate the electric fields in the 
piezoelectric layers and do not solve the coupled equations of piezoelectricity directly. Alternatively, 
(Tzou and Ye, 1996) proposed the analysis of piezoelectric shells as a layerwise assembly of 
curvilinear solid piezoelectric triangular elements, while Heyliger et al. (1996) developed discrete- 
layer mechanics and a finite element for laminated piezoelectric solids. 

The work in this paper attempts to remedy some of the problems in previous works (uncoupled 
approaches, limitations in piezoelectric laminations, large problem sizes, and so forth) while 
maintaining some of their advantages (formulation in curvilinear coordinates, layerwise approach, 
computational efficiency). Consequently, this paper presents the theoretical foundations of a new 
coupled laminate theory and, more importantlly, the development and evaluation of a corresponding 
finite element which enables the formal anlaysis of piezoelectric composite shells. The mechanics are 
truly coupled, that is, they involve approximate through-the-thickness fields for both displacements 
and electric potential, and are formulated using the coupled equations of piezoelectricity in the 
curvilinear coordinates. Moreover, the proposed laminate shell theory utilizes different types of 
approximations for the displacement and electric potential, that is, first-order shear theory type of 
assumptions for the displacements and the so-called discrete-layer (or layerwise) approximation for 
the electric potential (Heyliger et. al. (1994); Saravanos and Heyliger (1995, 1996)). The combination 
of mixed through-the-thickness approximations for the displacement and electric potential is a unique 
feature of the mechanics (termed, thereafter, as "mixed-field laminate shell theory") which enables 
the analysis of thin and moderately thick piezoelectric shells of general laminations with reasonable 
computational efficiency, while maintaining sufficient detail in the approximation of the electrical 
fields. Structural mechanics are also formulated, based on the laminate shell theory for the numerical 
analysis of piezoelectric composite shells and an 8-node curvilinear shell element is developed. 
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Numerical results for laminated composite plates are compared with results from an exact solution. 
Additional cases investigate the quasi-static and free-vibration response of an adaptive composite ring 
with a piezoceramic layer, and the quasi-static response of a [0/90/45/-45], semi-circular composite 
shell with discrete cylindrical piezoelectric actuators and sensors. 


2. Piezoelectric Laminated Shells 

This section describes the analytical formulation for curvilinear laminates with embedded sensory and 
active piezoelectric layers. The curvilinear laminate configuration is shown schematically in Fig. 1. 
Each ply of the laminate remains parallel to a reference curvilinear surface A*. An orthogonal 
curvilinear coordinate system O^T\^ is defined, such that the axes ^ and x\ lie on the curvilinear 
reference surface A^ while the axis £ remains straight and perpendicular to the layers of the laminate. 
A global cartesian coordinate system Oxyz is used to define A„ , hence, a point r=(x,y,z) on the 
curvilinear laminate is, 

C (1) 

where, T 0 =(x 0) y 0) z^ are the Cartesian coordinates of the reference surface A,,, and £ overhat indicates 
the unit vector perpendicular to the reference surface. 

Governing Material Equations 

Each ply is assumed to consist of a linear piezoelectric material with properties defined on the 
orthogonal curvilinear system O^ri^, and constitutive equations of the following form, 

J ' f 1 * d t (2 > 

D i m d b a j * z aPk 


or equivalently. 
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( 3 ) 


a i " C tf S J ~ e uFk 
Di - e£j ♦ 

where: ij = and k,l= 1,...,3; and S £ are the mechanical stresses and engineering strains in 
vectorial notation; is the electric field vector; is the electric displacement vector; s 5 and C fi are 
the elastic compliance and stiflness tensors; dy and e,j are different forms of the piezoelectric tensor; 
and % is the electric permittivity tensor of the material. Superscripts E, a, and S indicate constant 
electric field, stress and strain conditions, respectively. The axes 1, 2, and 3 of the material are parallel 
to the curvilinear axes r|, and £, respectively. The materials are assumed to be monoclinic class 2 
crystals with a diad axis parallel to the £ axis. The assumed material class is general enough, such that 
eqs. (1) or (2) may encompass the behavior of off-axis homogenized piezoelectric plies, as well as, 
passive composite plies. 


The tensorial strain S g and electric field components in a curvilinear coordinate system are related to 
the displacement and electric potential respectively as follows (Fung, 1965) 




(4) 


E C 


i a<j> 

Jzti dXi 


(5) 


where: gy is the metric tensor and are the Christofel Symbols of the curvilinear system; the 
overbar was used to indicate tensorial components in the curvilinear system, and implied in the rest 
of the paper. For the case of the curvilinear system 0£r|£ defined in Fig. 1, the Euclidean metric 
tensor takes the form (Soedel, 1993), 

e n -(Ui/Rfyn, * 33-1 ( 6 ) 

where R* are the local radi of curvature, and g° n , g 0 ^ are the components of the metric tensor on the 
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surface A,, (£=0) defined as 


Zn 


r 2 2 2 
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o I 2 2 2~ 
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(7) 


Consequently, after substituting in eq. (4) and with substantial derivations, the strain-displacement 
relationships become (see also, Soedel (1993)), 
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( 8 ) 


where u, v, w are displacements in the curvilinear system. Similarly, eq. (5) relates the electric field 
vector E* to the electric potential 4> as follows (similar to Tzou and Zhong (1993)), 


E x 4> 5 

E * —r 1 — +.n 

g^l.C« 2 ) 

E 3’ -+.C 


(9) 


Mixed Field Laminate Theory 

A new laminate theory is proposed in this paper which assumes linear displacement fields through the 
thickness of the laminate for the displacements u, v (displacement along the x and h axes respectively) 
while the transverse displacement w remains constant through the thickness (Fig. 2). However, a 
layerwise electric potential field is assumed through the laminate, consisting of N discrete continuous 
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segments (Fig. 2). Previous work by Saravanos and Heyliger on layerwise theories for piezoelectric 
plates has demonstrated the advantages and necessity of layerwise approaches in capturing the 
complicated electric fields and interactions which are present in piezoelectric actuators and sensors. 
Consequently, by mixing a layerwise electric potential field with first-order shear theory assumptions 
for the displacements, a new piezoelectric shell laminate theory is postulated (termed accordingly as 
"mixed-field laminate theory") which can: (1) capture the electric heterogeneity through-the- 
thickness induced by the embedded piezoelectric sensors & actuators, and (2) model thin and/or 
moderately thick piezoelectric shell laminates with minimal computational cost. 

The laminate is subdivided in N-1 piezoelectric sublaminates (or discrete-layers), and the subdivision 
can be arbitrarily controlled according to the configuration of the piezoelectric layers and/or the 
required detail of approximation. A continuous electric potential is assumed in each sublaminate, such 
that a C„ continuous variation results through the thickness of the laminate (see Fig. 2). The 
displacements and electric potential of the mixed-field theory take the following form, 

«(£,n»C,0- - CP g (M>0 

♦ CP 

(1®) 

■KC.'n.C.O- E ‘fr'CC.wMO 
>1 

where u°, v°, w° are displacements along the v\ and £ axes, respectively, on the reference surface 
superscript) indicates the points £ J at the beginning and end of each discrete layer; is the electric 

potential at each point £ (see Fig. 2); Y j (z) are interpolation functions; and ^ are the rotation 
angles defined as. 
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Two unique advantages of the laminate theory are obvious: (1) the complete electromechanical state 
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of the smart laminate is represented; and (2) the formulation entails the inherent option to select the 
detail of approximation of the electric field. Linear interpolation functions Y(z) were considered in 
this paper. 

The Love's assumption is further implemented, that is, the local radii of the shell are substantially 
higher than the thickness (h/R* « 1), yielding (1+ £/Rj * 1). The Love's assumption is appropriate 
for shallow shells, yet, additional work will be presented in the near future for deep shells with the 
Love's assumption removed. Now, in the context of eqs. (10), the engineering strains become, 

S 3 (M> C>0- o (12) 

Sfta.ifr i- 4,5 


where S° and k are the strain and curvatures at the reference surface, defined as follows: 
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( 13 ) 
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( 14 ) 


The electric field vector also becomes. 


9 




( 15 ) 


( 16 ) 


Equations of Motion 

The variational form of the equations of motion in the orthogonal curvilinear system is, 

f { -(6 H(s#y6T maiW C ♦ fibufr f**5ydr - o (17) 

where 8H is the variation of the electric enthalphy and ST the variation of the kinetic energy, defined 
as, 

6H - 6S i a l -6Ep i - bSpp f bE# t p f bS fl F r bE k * 1 p l 
6T - -6u i pjti i 

t and D overbar are respectively the surface tractions and electric displacement on the boundary 
surface F, 

D-Djiij *V-1,...,3 (19) 

[J] is the jacobian matrix of the transformation between the global Cartesian and curvilinear system. 
For the curvilinear system the determinant of the Jacobian takes the form: 

( 20 ) 

It is now possible to rearrange the variational statement and separate the through-the-thickness 
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integration, as follows: 

L { /o*‘C 6H ( 5,E >~ 6r >* + i+t>4>D)jr-0 (21) 

* 

The electric enthalpy and kinetic energy of the laminate may now be defined as, 

<6H l , 6T l >- f* <bff s bT> g\\g-n (22) 

where h is the laminate thickness. By combining eqs. (18,22) and integrating through-the-thickness, 
the variation of the electric enthalpy of the piezoelectric laminate is obtained as a quadratic expression 
of the generalized strain/electric field and the generalized laminate stiffness, piezoelectric and 
dielectric permittivity matrices, 

6H l - - 

-£ ( 6 S ;EZE?.6E?E^;.6kE"Er.&E~E"k> -£ £ iE^E," (23) 

m-l MM 


In the above equation, [A], [B], and [D] are the stiffness matrices of the curvilinear laminate, 


Ay By D s > - g° lf :n £ ‘J-IX6 


Ay ■ ? 11?22 £ » *J a 4, 5 , 


(24) 


[E m ] overbar and overhat are the piezoelectric laminate matrices, 


< e’ e” > . g°g£ E f C %T"j(0<l,C>4{. .-3, j. 14,6 

w * u Jf. 


M 

» * /• 




i-lJ2, >4,5 


(25) 


and [G ,nn ] are the laminate matrices of electric permittivity. 
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( 26 ) 


G,” - Znzi E I - '"** T’(C)T*(Oi{ 

M 

G 3 T - £ C Cu e^O^C)^ 

M ■»£/ 

where L is the number of plies in the laminate. All remaining matrix terms not shown above are zero. 

Combining equations (18, 22) and integrating through-the-thickness, the kinetic energy of the 
laminate takes the form, 

j> >1,2 (27) 

where u°f= {u°y,w°} and Pi={^, P n ); p A , p B , p D are the generalized densities, expressing the mass, 
mass coupling and rotational inertia per unit area, respectively, of the laminate, 

<P A ,P B ,P D > - E f S, (28) 

U\ 


3. Piezoelectric Shell Structures 


The formulation of the governing equations in the orthogonal curvilinear system and the attained 
separation of the through-the-thickness integration in the equation of motion (21) enables the 
development of structural solutions by using approximations of the generalized electromechanical 
state (displacements, rotation angles and electric potential) on the reference surface A^, of the 
following type. 


M 

«/«.«!. 0- E“/»^'«.l). 

(-1 

M 

E P!»^'«.i)>, >U (29) 

w 

M 

(•i 


where superscript i indicates the reference surface displacement, rotation angle and generalized 
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electric potential components corresponding to the i-th in-plane interpolation function N^x^y). This 
formulation forms the basis for an number of approximate solutions. Combining eqs. (13, 14, 16, 29), 

the generalized laminate strains and electric field at a point (£,ri) of the reference surface take the 

* 

form: 
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(32) 


Finite Element Formulation 


For structural problems with general boundary, geometry and material configurations, local 
interpolation functions may be used in eqs. (29) to develop finite-element based solutions. 
Substituting these approximations into the generalized equation of motion (21) and collecting the 
coefficients, as mandated by eqs. (23, 27), the governing dynamic equations of the structure are 
expressed in a discrete matrix form as. 


KJ o f {C7> 1 [*J [K u +1 f {U} 1 | {F(/» 

o oJ[ {$> j [[**] [r^Ji {4»> J"1 {2('» 


(33) 
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The above submatrices K* and are the elastic, piezoelectric and permitivity matrices and 
is the mass matrix. They are calculated from the generalized laminate matrices defined in eqs. (24-26) 
as determined by the variational statement. At each finite element they take the form: 


[* J* ■ f A 


m pi 

IS] [D]J 


[R *] </e*i 


trjjf' . I [r «f 


L 

WJ? - / [Jf'f 


[i~] 

|p”) 


[R 1 *] dldx\ 


[K^ - j [J? ^] r [G mt ]R * dldi\ 

[diag( p A )\ [diag{p B j\ 
[diag(p B y\ [duxgCp 0 )] 


[N*] dldxi 


(34) 


where: R" and R ra are the strain and electric field interpolation matrices defined by eqs. (30-3 1) and 
(32), respectively; N“ are the displacement interpolation functions defined by eq. (29). 


Based on the above formulation, an 8-node (M=8) finite element was developed with bi-quadratic 
shape functions of the serendipity family. The same interpolation functions are used to define the 
reference surface with respect to the coordinates of each node at the Cartesian system Oxyz 

\x V V N '(^> (* ** ** *> (35) 

M 


The radii R l3 R 2 and the coefficients g° n , g 0 ^ (eq. 7) of the reference surface are calculated in each 
element based on the previous parametric surface representation. 

Assuming that both sensory and active piezoelectric layers are embedded into the structure, the 
electric potential vector is subdivided in a free or sensory component d> F representing the voltage 
output at the sensors, and a forced or active component <j> A representing the voltage imposed on the 
active piezoelectric layers, such that {cj>} = { 4> F ; <& A }. Separating the active and sensory potential 
components in eq. (33), the equations take the following form 
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( 36 ) 


J o 

<“■> l . 

[*J 


1 lu) L 

0 °. 

l{* F il 



1 <* f >j 


{Q F (t)}-[K%U* A } 


where superscripts F and A indicate the partitioned submatrices in accordance with the selected 
sensory and active configuration, respectively. The left-hand side includes the unknown 
electromechanical response of the structure {u, <|> F }, that is, the resultant displacements and voltage 
at the sensors. The right-hand includes the excitation of the structure in terms of mechanical loads 
and applied voltages on the actuators. The electric charge at the sensors Q F (t) remains constant with 
time (practically open-circuit conditions) and is assumed equal to zero. 


Among the obvious advantages is the capability of the mechanics to model the response of the 
piezoelectric structure either: in ’'active" mode, that is, with specified voltages Acj> A applied across 
the piezoelectric layers to induce a desirable deflection/strain state; or in "sensory" mode where 
displacements or mechanical loads are applied to the structure and the resultant voltage or charge is 
monitored; or in combined "active/sensory" mode. The above dynamic system may be solved to 
obtain, either the modal characteristics (free-vibration), or the forced frequency response, or the 
transient response of the piezoelectric shell. 


4. Evaluations and Discussion 

Applications and evaluations of the developed mechanics on various piezoelectric structures are 
presented in this section. As in previous publications, the standard laminate notation is expanded 
such that piezoelectric layers are indicated with the letter p. Results are presented for simply 
supported [p/0/90/0/p] cross-ply Graphite/Epoxy composite plates with attached piezoceramic 
(PZT-4) layers. Results are also presented for two cylindrical shell structures: (1) a closed 
cylindrical titanium shell [Ti/p] with a continuous piezoceramic (PZT-4) layer on the outer 
surface; and (2) a cantilever semi-circular laminated [0/90/45/-45], Graphite/Epoxy shell with 
distributed cylindrical piezoceramic (PZT-4) patches on the inner and outer surface. The 
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properties of the materials are provided in Table 1. 


[p/0/90/0/p] Hybrid Plate 

Results are presented for the ffee-vibration response of simply supported [p/0/90/0/p] cross-ply 
Graphite/Epoxy composite plates with surface attached piezoceramic (PZT-4) layers. The electric 
potential on the outer (free) surface of the piezoelectric layers was always forced to remain zero. 
The electric potential in the inner surface may be either free or forced to remain zero, resulting in 
two different electric boundary conditions, termed as open- and closed-circuit conditions 
respectively. The thickness of each composite ply and piezoelectric layer are 0.267h and 0. lh 
respectively, where h is the plate thickness. The composite plies and piezoelectric material were 
assumed to have equal mass density in this case. The response of these plates has been previously 
studied via exact solutions (Heyliger et. al. 1995) and discrete-layer mechanics (Saravanos and 
Heyliger, 1996), hence, their analysis provides excellent assessments of the accuracy and range of 
capabilities of the present mechanics . 

Table 2 shows predictions of the fundamental natural frequency for a moderately thin plate 
(a/h=50) and a thick plate (a/h=4) with either open or closed-circuit conditions imposed on the 
piezoelectric layers. The case of interest here is the moderately thin plate because the mechanics 
are developed for thin and moderatelly thick shells; whereas, the thick plate example is an extreme 
case presented to quantify any limitations of the present approach. In the case of the t hi n plate, the 
predicted natural frequencies converge, either as the in-plane mesh density is increasing (M) or the 
number of discrete layers for the electric potential through-the-thickness (N). More importantly, 
the predicted natural frequency agrees well with the exact piezoelasticity solution. As expected, 
lower frequency values were predicted for closed (C), than with open (O) circuit conditions. Figs. 
3a-b shows the predicted modal displacements and electric potential fields through-the-thickness 
of the thin plate, and as seen, excellent agreement was obtained between the present approach and 
the exact solution. Although the present method is not indended for thick shells, results for the 
case of a thick plate are also presented in Table 2 and Fig. 4. The predicted frequency values in 
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Table 2 overestimate the exact solution results by 12.7% and 21.5% for (O) and (C) conditions, 
respectively, as the laminate theory does not capture the layerwise variations in the displacement 
fields (Fig. 4a-b). However, Fig. 4c shows that the present method did capture the layerwise 
variations of the modal electric potential and yielded excellent agreement with the exact solution. 

Closed Cylindrical Shell 

Results for the quasi-static and free-vibration response of an active/sensory cylindrical ring are 
presented in this section. The ring consisted of one 3 mm thick Titanium layer and one 
continuous piezoceramic layer of PZT-4 (1 mm thick) attached on the outer surface. The inner 
radius was 289 mm, and the width was 348 mm. Fig. 5 shows one quarter of the ring, the laminate 
configuration and the curvilinear coordinate system which was used. Symmetry conditions were 
applied at two planes of symmetry (at 0= 0 and 90° ) in the circumferential direction (£) , and at 
the plane of symmetry in the axial direction (r\) in order to reduce the size of the model. Two 
discrete layers (N=3) were used for the electric potential field through-the-thickness (0, one f° r 
the piezoelectric and one for the titanium layer. The electric potential in the inner surface of the 
piezoelectric layer was always equal to 0 Volts. Two types of loading conditions were examined: 
(1) application of a sinusoidal electric potential on the outer surface of the piezoelectric layer; (2) 
application of mechanical load with free electric potential on the outer piezoelectric surface. 
Apparently, the previous two loading conditions were selected to investigate the active and 
sensory response, respectively, of the cylindrical shell. 

Active Ring. The radial displacements at 0= 0 and 90° at the center of the ring (along the r\ 
axis), induced by a sinusoidal electric potential (1OOcos20 Volts) applied on the outer surface of 
the piezoelectric actuator, are shown in Table 3 for various mesh densities. The first five natural 
frequencies of the active ring are also shown in Table 3 which correspond to closed-circuit 
electric conditions. The induced displacements and the natural frequencies seem to converge well 
as the mesh density increases. There is also very good agreement between the present 
predictions and numerical results obtained with a recently developed layerwise finite-element for 
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laminated piezoelectric solids (Heyliger et. al, 1996), which validates the accuracy of the present 
method. Fig. 6 shows the induced active radial deflections along the circumferential direction at 
the center and the free edge of the active ring. The applied sinusoidal electric potential produces 
an almost pure ovalization of the active shell. Considering that the applied voltage in this example 
was approximately 4 times lower that the maximum allowed voltage on the actuator, the 
maximum active deflection which may be effectively induced is not insignificant. Such ovalization 
capabilities may be important to many aeronautical applications. 

Sensory Ring. This example examines the behavior of the piezoelectric layers as a distributed 
sensor. The electric potential on the outer surface of the piezoelectric layer remains free, while a 
radial line load F^= 656 N/m was applied at 0= 90° . As in the previous case, induced deflections, 
sensory potential and natural frequencies of the sensory ring are shown in Table 4, together with 
convergence studies and correlations with other published results. Again the convergence of the 
present study is excellent, with the sensory electric potential requiring the higher mesh density. 

The agreement between the present method and layerwise piezoelectric solid mechanics (Heyliger 
et al., 1996) ranges from excellent (natural frequencies) to fair (displacements). The variations of 
the radial deflection and corresponding sensory electric potential in the hoop direction are shown 
in Fig. 7. The application of the line load also results in the ovalization of the ring, and the sensory 
potential has a nearly sinusoidal form. The results in Fig. 7 demonstrate the usefulness of the 
present mechanics in predicting the not so apparent sensory voltage patterns in distributed 
piezoelectric sensors which correspond to various deformed shapes of a smart shell structure. The 
usefulness of these predictions in estimating the minimum number and location of discrete 
piezoelectric patches is also apparent. Yet, as shown in the next case, the method has the 
capability to explicitly analyze the response of active/sensory composite shells with discrete 
piezoelectric actuators and sensors. 

[0/90/+45/-45], Cantilever Semi-Circular Shell 

The response of a semi-circular cantilever [0/90/+45/-45], Graphite/Epoxy smart shell shown in 
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Fig. 8 was also analyzed. Four curved piezoceramic (PZT-4) patches are attached on each side of 
the composite shell, covering 80% of the free area. The mid-surface radius is 291 mm and the 
width of the shell (L^) is 152.4 mm. The thickness of each composite ply is 0.12 mm, while the 
thickness of each piezoelectric patch is 0.24 mm. The curvilinear coordinate system and the finite 
element discretization are also shown in Fig. 8b. Three discrete layers (N=4) were used for the 
electric potential field through-the-thickness (Q, one for each piezoelectric layer and one for the 
composite laminate. In all cases, the electric potential at the bonded surface of the piezoelectric 
patches was forced to remain zero. The specification of the other two generalized electric 
potential values can result in various active and/or sensory configurations of the smart shell. Two 
such active and sensory configurations were analyzed. 

Active Case. The induced active deflections in the radial (w°) and hoop direction (u°) are shown 
in Fig. 9. All piezoceramic patches were configured as actuators, with 100 Volts applied on their 
free surfaces, thus resulting in electric fields of opposite polarities at the piezoceramics of the 
inner and outer surface. For comparison purposes, the radial and circumferencial deflections of 
the shell with continuous piezoceramic layers of same thickness are also shown in Fig. 9. 
Although, the calculated deflections with either the discrete or the continuous piezoceramic 
actuators seem to be in good agreement, the effect of discrete actuators on the deflected shape of 
the shell is obvious. These effects can be attributed to the non-uniform actuation and the non- 
uniform laminate stiffness resulting from the piezoelectric patches. It is important to point out, 
however, the relatively high radial and hoop deflections at the free-end (approximately 5 mm) 
obtained in this case. It highlights the possibility to achieve substantialy large, yet, accurate and 
rapid free-end positioning with such active cantilever rings. 

Sensory Case. The response of the cantilever shell with all piezoelectric patches configured as 
sensors (free electric potential) was also modelled. The more realistic case of continuous 
electrodes on the free surface of each sensor was also modeled. A continuous surface electrode 
forces a uniform electric potential over the surface, which was represented with equality 
constraints on the electric potential. A line load in the hoop direction F^ = 656 N/mm was applied 
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at the free end. The resultant displacements and corresponding voltages at each sensor are shown 
in Figs. lOa-b, respectively. The predicted response in the case of continuous piezoelectric layers 
is also shown. The presence of discrete piezoelectric sensors has a definite effect on the deflected 
shape. It is believed that the effects on the deflected shapes have resulted in the more dramatic 
differences between the voltage of discrete vs. continuous sensors. Fig. 10 also illustrates that the 
response of discrete sensors may be only roughly approximated with the consideration of a 
continuous layer. Hence, many proposed methods which attempt to replicate discrete sensors with 
segments of a continuous piezoelectric layer may lead to considerable error. 

5. Summary 

Mechanics for the analysis of laminated composite shells with piezoelectric actuators and sensors 
were presented. The mechanics were based on a new mixed-field theory for curvilinear laminates 
which combines single-layer assumptions for the displacements with a layerwise representation of 
the electric potential. The resulting governing coupled equations for piezoelectric shell laminates 
and structures were developed in curvilinear coordinates. Based on them, finite-element based 
analysis procedures for piezoelectric composite shell structures were described, and an 8-node 
shell finite-element was formulated for the static and dynamic analysis of composite shells 
containing piezoelectric actuators and sensors. The described mechanics, the finite element and 
the computational procedure were encoded in prototype software. 

Excellent correlations with previously reported results for thin piezoelectric plates were obtained. 
Additional evaluations and comparisons were reported for the active and sensory response of a 
titanium cylindrical ring with a continuous piezoceramic actuator or sensor, and the capability to 
actively ovalize closed cylindrical shells was illustrated. Finally, the response of cantilever semi- 
circular laminated Graphite/Epoxy shells with discrete piezoelectric patches was also analyzed. 
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Table 1 Mechanical Properties 
(e 0 =8.85 1 O ' 12 farad/m, electric permittivity of air) 


Property 

Titanium 

Gr/Epoxy 

PZT-4 

Elastic Properties: 




E u (GPa) 

114 

132.4 

81.3 

Eb (GPa) 

114 

10.8 

81.3 

E 33 (GPa) 

114 

10.8 

64.5 

G* (GPa) 

43.8 

3.6 

25.6 

G u (GPa) 

43.8 

5.6 

25.6 

G 12 (GPa) 

43.8 

5.6 

30.6 

Vl2 

0.3 

0.24 

0.33 

V i3 

0.3 

0.24 

0.43 

Va 

0.3 

0.49 

0.43 

Piezoelectric 
coefficients (l(T 12 m/V): 




d 3i 

0 

0 

-122 

d 32 

0 

0 

- 122 . 

d 2 4 

0 

0 

495. 

d 15 

0 

0 

495. 

Electric Permittivity: 




e n /e„ 

1475. 

3.5 

1475. 


1475. 

3.0 

1475. 

We. 

1475. 

3.0 

1300. 

Mass Density p (kg/m 3 ) 

2768. 

1578. 

7600. 
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Table 2. Predicted Natural Frequencies of [p/0/90/0/p] square simply supported plate. 
(C)- closed circuit; (O)- open circuit 


Mesh 

Density 

Discrete 

Layers 

(N-l) 

Fundamental Frequency f t (F/h p/ n 
Iff Hz (kg/m) 1/2 



oi/h~50 (C) 

a/h=50 (O) 

a/h=4 (C) 

a/h = 4 (O) 







4x4 

3 

236.390 

258.673 

255.052 

176.576 

8x8 

3 

232.574 

254.969 

163.594 

176.516 

12x12 

3 

231.938 

254.379 

163.593 

176.609 







8x8 

7 

232.647 

255.036 

163.849 

176.713 

8x8 

20 

232.665 

255.052 

163.909 

176.760 







Exact: 


245.941 

245.942 

145.339 

145.377 
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Table 3. Predicted displacements and natural frequencies of Active [Ti/PZT-4] 
cylindrical ring with 1OOcos20 Volts applied on the piezoelectric layer. 



Mesh Density 

0 

Heyliger et 
al. (1996) 

( 1 6x2,N=5) 


5x1 

10x1 

15x1 

20x1 

15x2 

15x4 

Radial 
Displacement 
at Center 
(mm) 








w(0=O) 

0.0318 

0.0976 

0.1103 

0.1130 

0.1084 

0.1065 

0.1112 

w(0=9O) 

-0.0318 

-0.0976 

-0.1103 

-0.1130 

-0.1084 

-0.1065 

-0.1112 









Natural 

Frequencies 

(Hz) 








fi 

57.98 

33.01 

31.06 

30.68 

30.96 

30.88 

31.27 


258.39 

174.00 

167.64 

166.27 

167.11 

166.68 

170.42 

f 3 

640.39 

418.22 

399.73 

395.23 

398.67 

397.70 

407.29 


1250.39 

776.09 

730.57 

718.47 

728.73 

726.98 

745.21 

f 5 

2093.49 

1263.63 

1166.20 

1139.05 

1163.49 

1160.68 

1190.48 
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Table 4. Predicted displacements and natural frequencies of Sensory [Ti/PZT-4] 
cylindrical ring with a radial line force F^= 656 N/m applied at 0= 90° . 



Mesh Density 


* 

Heyliger et 
al. (1996) 
(16x2,N=5) 


5x1 

10x1 

15x1 

20x1 

15x2 

15x4 

Radial 

Displacement 
at Center 
(mm) 








w(6=0) 

-1.004 

-2.960 

-3.312 

-3.387 

-3.348 

-3.387 

-3.750 

w(0=9O) 

1.128 

3.233 

3.607 

3.687 

3.647 

3.687 

4.083 









Electric 
Potential at 
Center (V) 








4> F (0=O) 

68.90 

231.80 

266.60 

274.60 

285.30 

284.40 

N/A 

<f(e=90) 

-107.50 

-361.30 

-428.80 

-453.10 

-458.30 

-457.70 

N/A 









Natural 

Frequencies 

(Hz) 








f, 

58.83 

34.51 

32.66 

32.28 

32.48 

32.31 

32.60 

f 2 

264.06 

182.15 

176.07 

174.76 

175.19 

174.35 

175.11 

f 3 

654.90 

437.21 

419.13 

414.72 

417.29 

415.35 

417.92 


1278.99 

811.45 

766.07 

754.85 

762.27 

758.54 

763.94 

fs 

2132.06 

1320.94 

1223.55 

1196.25 

1217.60 

1211.25 

1220.01 
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Figure 1 . — Curvilinear piezoelectric laminate and coordinate systems. 
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(a) Smart Laminate (b) Piezoelectric Shell Laminate theory assumptions 


Figure 2. — Typical piezoelectric laminate configuration . (a) Concept ; (b) Assumed through-the-thickness 
displacement and electric potential fields. 
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Normalized Thickness (z/h) 



Figure 3. — Through-the-thickness modal displacement and electric potential distributions for a thin [p/0/90/0/p] 
plate (a/h=50). Fundamental mode. 
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Figure 5. — Quarter model of the cylindrical ring with an attached piezoceramic layer. 


Radial Deflection, (m) 


Ti/PZT-4 Ring 

100 cos2t Volts applied 
15x2 mesh, 1/4 model 



Center 


Free-edge 


0 15 30 45 60 75 90 


Hoop angle (degs) 


Figure 6. — Actively induced radial deflection of the cylindrical ring with 100cos29 Volts sinusoidal electric 
potential applied on the piezoelectric actuator. 
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Center 

Free-edge 


0 15 30 45 60 75 90 


Hoop angle (degs) 



Figure 7. — Radial deflection and corresponding sensory electric potential of the cylindrical ring with a radial line 
load of 656 N/m applied at 0=90°. 
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Deflection (m) 
(Times IE-3 ) 


* 



Figure 9. — Active radial and circumferencial deflections of the semi-circular cantilever shell with 100 Volts 
applied on each piezoelectric patch. 
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